Quadrupole collective modes in trapped finite-temperature Bose-Einstein condensates 
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Finite temperature simulations are used to study 
quadrupole excitations of a trapped Bose-Einstein conden- 
sate. We focus specifically on the m = mode, where a 
long-standing theoretical problem has been to account for an 
anomalous variation of the mode frequency with temperature. 
We explain this behavior in terms of the excitation of two sep- 
arate modes, corresponding to coupled motion of the conden- 
sate and thermal cloud. The relative amplitudes of the modes 
depends sensitively on the temperature and on the frequency 
of the harmonic drive used to excite them. Good agreement 
with experiment is found for appropriate drive frequencies. 

PACS numbers: 03.75.Fi, 05.30. Jp, 67.40. Db 



One of the most significant challenges facing theorists 
studying Bose-Einstein condensation (BEG) is to under- 
stand the finite-temperature properties of trapped gases, 
where the noncondensed fraction is significant. For in- 
stance, while solutions of the Gross-Pitaevskii (GP) equa- 
tion for low-lying collective modes yield good agree- 
ment with experiments at low temperatures the be- 
havior at elevated temperatures is still not fully under- 
stood. The quadrupole modes studied in an early exper- 
iment by Jin et al. have been especially problematic. 
In this experiment, a harmonic modulation of the trap 
frequency was used to selectively excite either m = or 
TO = 2 modes, and the frequency and damping of the 
modes as a function of temperature, T, were extracted 
from the subsequent response. Although several theories 
have been proposed that account for the observed damp- 
ing and frequency of the m ~ 2 mode [^-0, the 
experimental data for the to = frequency has yet to 
find a satisfactory explanation. 

The difficulty is apparently related to the fact that 
the frequency of the m = condensate mode lies rela- 
tively close to that of the corresponding noncondensate 
mode. In particular, the condensate frequency at T = 
in the Thomas- Fermi (high density) limit is cj = l.SOtJ^ 
[Ql (where u}± is the axial trap frequency), while the 
thermal cloud has a natural frequency of oscillation at 
u) — 2llJj_. It should be emphasized that the latter is not 
a true collective mode but is rather a coherent motion 
of the atoms in a harmonic potential. At temperatures 
below Tc, where the condensate coexists with a signifi- 
cant thermal fraction, these two modes of oscillation are 
strongly coupled by mean- field interactions. As a result, 
a harmonic drive will tend to excite a superposition of 
the two, and the observed m = condensate frequency 
will reveal this coupling. A consistent theoretical model 



should therefore include the full dynamics of the thermal 
cloud, and treatments that neglect fluctuations of the 
noncondensate |^ or include them perturbatively 
will not capture this behavior. 

To avoid this shortcoming, Refs. [ p^|jTl| ] used a vari- 
ational scheme to solve coupled GP and Boltzmann ki- 
netic equations to describe the dynamics of, respectively, 
the condensate order parameter and the noncondensed 
component. Normal modes corresponding to in-phase 
and out-of-phase oscillations of the two components were 
found, and the authors suggested that the m = exper- 
imental data could be explained as a cross-over between 
the two modes. As we shall see, this picture in fact has 
some validity, however a quantitative interpretation of 
the experiments requires a more accurate treatment of 
the thermal cloud. In particular. Landau damping, which 
is the dominant damping mechanism in this regime, was 
not included since the variational ansatz used precludes 
coupling to higher order modes. Moreover, in view of our 
earlier comment about the nature of the thermal cloud 
oscillation, it is unclear to what extent the simple normal 
mode picture for the thermal cloud is meaningful. 

In this work we model the experiment of Jin et al. 
Q using numerical simulations based upon the ZNG for- 
malism . This approach was previously applied to the 
scissors mode |l3| and yielded good agreement with ex- 
periment We evolve a GP equation simultaneously 
with a semiclassical kinetic equation, where the latter 
is solved by representing the thermal cloud by means 
of an ensemble of test particles which suffer collisions 
with each other as well as with the condensate. Impor- 
tantly, since the particles are also coupled to the conden- 
sate through mean- field interactions, our model includes 
Landau damping (and corresponding frequency shifts) of 
condensate oscillations. Further details of the numerical 



method can be found in |15|. 

By using different initial conditions we can separately 
excite motion of the condensate and thermal cloud. Al- 
ternatively, by imposing a harmonic modulation of the 
trap frequencies over a finite interval (in a similar fash- 
ion to experiment), we excite both components simulta- 
neously. At high temperatures we find that the conden- 
sate exhibits oscillations with two frequency components 
— one corresponding to a damped oscillation of the con- 
densate interacting with a quasistatic bath of thermal 
particles, and a second arising from the mechanical cou- 
pling to the oscillation of the thermal cloud. The relative 
amplitudes of the two components depends on temper- 
ature and on the frequency of the harmonic drive. By 
using different drives we can closely reproduce the ob- 
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served temperature dependence of the m = frequency. 

The GP equation for the condensate wavefunction, 
$(r, t), and the Boltzmann kinetic equation for the ther- 
mal cloud phase space density, /(p, r, t), are respectively 
given by 
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V/-VC/eff • Vp/ = Cl2[/]+C22[/]. (2) 



The above equations are derived under certain approx- 
imations jl2|. Importantly, thermal excitations are 
treated within the Hartree-Fock (HF) and semiclassi- 
cal approximations [ p6| , where they can be identified 
with particles moving in an effective potential Ucs — 
mLo']_{x^ + 2/^ + \^z'^)/2 + 2g{nc + n). Here, ^^(r, t) = 
|$(r, i)p and n{r,t) = /(dp//i^) /(p, r, i) are the con- 
densate and noncondensate densities respectively {g — 
Anh'^a/'m, where a is the s-wave scattering length and 
TO is the atomic mass). The integral C12 represents col- 
lisions that transfer atoms between the condensate and 
noncondensate, while C22 collisions take place between 
thermal particles only. The former process couples to the 
GP equation through R{r,t) = {h/2nc) /(dp/Zi^) Ci2[/] 
and leads to a change in the number of condensate atoms. 

Our simulation parameters are chosen to match those 
of the experiment where wj_ = 27r x 129Hz and 
A = Vs. The equilibrium densities of the condensate 
and noncondensate arc found using a self-consistent semi- 
classical procedure at each temperature. Due to 
evaporative cooling the total number of atoms, N, in the 
experiment varied with temperature. This is accounted 
for in our simulations by using experimental data for the 
temperature and number of condensate atoms, the vari- 
ables believed to be the most precisely determined [^s). 
Alternatively, we have also used the temperature and to- 
tal number of atoms to specify the state of the system. 
Although this leads to slightly different equilibrium con- 
ditions, it does not affect the general conclusions of our 
work. These experimental uncertainties, however, should 
be kept in mind when making comparisons between the- 
ory and experiment. 

We can excite quadrupole modes by imposing velocity 
fields of the form Vc = Aco{x,y) + Ac2{x,—y) onto the 
condensate, and v„ = AnQ{x,y) + An2ix,—y) onto the 
noncondensate. In the case of the condensate, the veloc- 
ity is established by multiplying the equilibrium conden- 
sate wavefunction by an appropriate phase factor. For 
the thermal cloud, we simply add the position-dependent 
v„ to the initial velocity of each test particle. Since the 
TO = and TO = 2 modes are decoupled from each other 
due to their different symmetries, it is advantageous to 
excite both simultaneously (e.g., Aco,^c2 > 0). The 
corresponding modes are then projected out in the en- 
suing evolution by evaluating the moments Q{n,c}o 



(to = 2). We note that the (Z = 0, to = 0) monopole 
and (Z = 2, TO = 0) quadrupole modes of an isotropic 
trap are coupled by the trap anisotropy [Q, giving rise 
to low and high frequency m — modes. Only the low- 
lying TO = out-of-phase breathing mode is appreciably 
excited by the initial conditions considered here. 

For our first simulations, we excite oscillations in the 
condensate alone (i.e., Ano = A^i = 0). The Qn 
amplitudes remain small during the course of the evo- 
lution, while the condensate oscillations are damped. 
To extract frequencies and damping rates we fit a sin- 
gle exponentially-decaying sinusoid to the data over the 
timescale of the simulation, uj\t = 25. The open 
squares in Fig. |l| show the frequencies as a function 
of the reduced temperature T' = T/T°, where T° = 
0.94ftwj_(AA^)^/"^/fcB is the ideal gas critical temperature. 
Apart from the highest temperature at T' — 0.9, the sin- 
gle function fit is very good, consistent with an excita- 
tion of just one mode of each symmetry. For both the 
m — and m — 2 modes we find a marked downward fre- 
quency shift with increasing temperature, in agreement 
with previous theories that neglect the full nonconden- 
sate dynamics l^-Q] , and with the out-of-phase "conden- 
sate" mode in Ref. Also, as in earlier treatments, 
we find that there is much better agreement with the ex- 
perimental TO = 2 data than with the to = data above 
T' = 0.6. 
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FIG. 1. Frequencies of the (a) m = and (b) m — 2 modes 
as a function of reduced temperature T' = T/T^. Data from 
the experiment of Ref. |^ is plotted with error bars, while our 
frequency results are shown for the condensate (open symbols) 
and noncondensate (filled symbols). Different initial condi- 
tions are represented by squares {A„o, An2 ~ 0) and triangles 
{A„o,A„2 > 0), where for both Aco,Ac2 > 0. 

We next investigate the influence of the thermal cloud 
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dynamics by exciting both components simultaneously. 
In these simulations, the velocity fields imposed on the 
condensate and thermal cloud are identical. The trian- 
gles in Fig. 1 show results of a single mode fit to the con- 
densate and noncondensate data. At high temperatures 
the thermal cloud oscillates with a frequency close to the 
ideal gas value, uj = 2uj_, which decreases as the tem- 
perature is lowered. The m = 2 condensate frequency is 
largely independent of the initial condition (indeed it is 
also unchanged for a harmonic modulation of the trap 
discussed below), apart from the highest temperature 
point where the condensate mode becomes ill-defined. In 
contrast, the m = condensate frequencies are gener- 
ally higher than those found by exciting the condensate 
alone, and approach the noncondensate value just below 
the transition. 

Further insight into this behavior can be gained by 
fitting two damped sinusoids to the m = condensate 
data. This shows that two modes are present at high 
temperatures, with one frequency close to that found for 
the condensate-only excitation. The second matches the 
noncondensate value, and indicates that the condensate 
is being driven by the thermal cloud at its own natu- 
ral frequency (through mean field interactions and C12 
collisions). Following Ref. [0 we refer to these as the 
"condensate" and "thermal cloud" modes, respectively. 
The observed condensate oscillation is thus a superpos- 
tion of the two modes, with the relative amplitudes be- 
ing sensitive to the initial conditions. By using a single 
mode fit to analyze the data, a frequency intermediate 
between the two mode frequencies is obtained, as shown 
by the open triangles in Fig. 1(a). In the case of the 
m = 2 modes, the condensate and noncondensate fre- 
quencies are sufficiently separated for the noncondensate 
to have a minimal effect on the condensate, and the lat- 
ter thus oscillates at a frequency which is independent of 
the method of excitation. 

Qualitatively, the behavior of the m — mode shown 
by the open triangles in Fig. hj(a) is reminiscent of 
that exhibited experimentally ff. However, a quan- 
titative comparison to experiment requires simulations 
which more faithfully reproduce the harmonic excitation 
scheme employed emperimentally. In the simulations to 
be described next, the axial trap frequency is modulated 
according to tf^(i) ~ uj'j_{l + €smUt) over a time interval 
of ui±t = 30, before allowing the sample to evolve in the 
original trap potential for another uj±t = 30. This ex- 
cites m = type oscillations in both the condensate and 
thermal cloud. Typical time-dependent plots of Qco and 
Qno are shown in Fig. || for = 1.95 and T' = 0.8. In 
the case of the condensate, we clearly see beating between 
the condensate and thermal cloud modes after removal of 
the excitation. This should be observable if experimental 
measurements are extended to longer times. 

In the experiments, a single mode fit was used to ex- 
tract mode frequencies from the data in a time interval 
of approximately uj±t = 15 following the excitation. We 
have followed this procedure by analyzing our data within 



the observation window indicated by the vertical lines in 
Fig. 1^. Very different results for the condensate frequency 
are obtained depending on the drive frequency which 
affects the relative amplitude of the condensate and ther- 
mal cloud oscillations. Fig. |^ summarizes our results by 
showing frequencies from single-mode fits at each T for 
drive frequencies in the range [1.75, 2.00] uj±. At low T 
the noncondensed component is small and has a minimal 
effect. As a result, all the frequencies are close to the 
condensate mode results shown by the open squares in 
Fig. ^ In contrast, above T' ~ 0.6 there is a significant 
spread in the frequencies extracted from the fits. Taking 
T' = 0.8 as an example, the condensate mode is excited 
when Q ~ 1.75luj_, while the frequency is close to that of 
the thermal cloud mode when f2 ~ 2ll!±. 
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FIG. 2. Plot of the condensate and noncondensate m = 
quadrupole moments QcO, QnO (in arbitrary units) as a func- 
tion of time, at T' = 0.8. During the first u)j_t — 30 (up to the 
first vertical line) the system is driven by a trap modulation 
of frequency Q = 1.95t<;_L and amplitude e — 0.02. The subse- 
quent evolution takes place in a static trap, where the second 
vertical line indicates the approximate range of experimental 
measurements. 

It is clear from our data that there are two distinct 
branches above T' ~ 0.7. The lower branch follows 
closely the condensate mode in Fig. |l| while the upper 
branch moves up to the thermal cloud frequency. Al- 
though the precise drive frequency used in the experi- 
ments is no longer known [ p8[ , the frequency was chosen 
to maximize the amplitude of the condensate oscillation. 
Using this criterion, we have identified the drive frequen- 
cies that give a large response by the solid points in Fig. ||. 
Again there are two clear branches, and although ex- 
periment could feasibly have followed either, the upper 
branch clearly seems most relevant. Given the uncer- 
tainties in the experimental conditions the agreement is 
pleasing, and is strong evidence that a cross-over from 
one branch to the other is responsible for the observed 
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temperature dependence. Our simulations nevertheless 
indicate that the lower branch should be observable if 
driven at the appropriate frequency. 
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FIG. 3. (a) Frequency of the m — mode as a function 
of temperature, with the experimental results of Q plotted 
with error bars. We display simulation results for six drive 
frequencies, with Q — 1.75tJx (circles), Q — 1.80a;x (squares), 
il = 1.85tJj_ (triangles), Q = 1.90tJx (stars), SI — 1.95ajx 
(diamonds), and fl = 2.00tJx (inverted triangles). The filled 
symbols at each T' represent the drives that produce the two 
largest condensate responses, (b) Damping rate versus T' , 
where our results (open symbols) are compared to experiment 
(closed symbols) , for the m = (triangles) and m — 2 (circles) 
modes. 

For completeness, we end by comparing our results 
for the damping rate against experiment. This is per- 
haps less instructive, since the experimental data con- 
tains large error bars. There are also uncertainties in the 
damping we obtain at the higher temperatures due to 
sensitivity of the results to both the form of the drive 
and the timescale of the fit. However, the scatter in 
the results at low temperatures is much smaller and, as 
Fig. ^(b) shows, the damping rates obtained in this tem- 
perature range are in good agreement with experiment. 

In summary, we study quadrupole oscillations of a 
Bose condensed gas at finite temperatures using simu- 
lations based on the ZNG formalism |l2j. This allows 
us to observe the coupled dynamics of the condensed 
and noncondensed components of the gas, which are par- 
icularly important when considering the m = mode 
where the characteristic frequencies of the two compo- 



nents are close. Simultaneous excitation of the conden- 
sate and thermal cloud by a harmonic modulation of the 
trap leads to a cross-over between the two frequencies at 
high temperatures, in accordance with experimental ob- 
servations. This agreement, taken with previous studies 
of the scissors mode indicates that a semiclassical 
Hartree-Fock description is adequate in describing finite- 
temperature dynamics in the coUisionless regime. 
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